Topological spectral correlations in 2-d disordered systems 



Vladimir E. Kravtsov 1 and Vladimir I. Yudson 2 
1 International Centre for Theoretical Physics, P.O. Box 586, 34100 Trieste, Italy 
and Landau Institute for Theoretical Physics, Russian Academy of Sciences, 2 Kosygina str., 117940, Moscow r-n, Russia 
2 Institute of Spectroscopy, Russian Academy of Sciences, 142092, Troitsk, Moscow r-n, Russia. 

(February 1, 2008) 



00 ■ 

o\ ■ 

o> : 

Oh! 
<D . 

m ■ 
t> ■ 



> 
in 
m 

(N 

O 

00 
ON 



c 

o 
o 



It is shown that the tail in the two-level spectral correla- 
tion function R(s) for 2D disordered systems depends on the 
system geometry and the boundary conditions. In particular, 
for closed surfaces (with no boundary) R(s) — — x/(67r 2 /3s 2 ), 
where /3 = 1, 2 or 4 for the orthogonal, unitary and symplec- 
tic ensembles, and x — 2(1 —p) is the Euler characteristics of 
the surface with p "handles" (holes). The result is valid for 
g <C s < g 2 for j3 = 1, 4 and for g -C s < ff 3 for ,3 = 2, where 
(? ^ 1 is the dimensionless conductance. 

PACS numbers: 73.20.Fz, 71.25.Mg, 72.15.Rn 



Spectral correlations in complex and disordered sys- 
tems are closely related with the basic symmetries of a 
system, the nature of its dynamics and the structure of 
the corresponding eigenstates. It is remarkable that there 
are few universal spectral statistics which all the generic 
complex and disordered systems obey in the thermody- 
namic (TD) limit. According to the well known Bohigas- 
Giannoni-Schmidt conjecture jjj, a generic chaotic sys- 
tem is described by the Wigner-Dyson (WD) spectral 
statistics which follows from one of the three universal 
Gaussian ensembles of random matrices [||. In contrast, 
all the generic integrable systems obey the Poisson spec- 
tral statistics. Depending on the strength of disorder, 
spectral statistics of a d-dimensional disordered system 
of non-interacting particles also flow to one of the above 
universal statistics as the system size L tends to infin- 
ity. For strong disorder and for low-dimensional systems 
(d = 1,2) where all states are localized, the spectral 
statistics in the TD limit are Poissonian. For weakly 
disordered 3D systems where all states are extended, the 
spectral statistics are identical to the WD || in the TD 
limit. 

However, the disordered systems with d > 2 are known 
to undergo the Anderson metal-insulator transition at a 
certain critical disorder. At this point the localization 
(correlation) length £ diverges and the spectral statistics 
should be independent of the system size L, provided that 
the energy difference u> = sA is measured in units of the 
mean level spacing A oc L~ d . Thus at the Anderson 
transition there exists a special fixed point, the critical 
spectral statistics (CSS) 

An amazing property of some CSS found recently Q is 
the sensitivity to boundary conditions and a shape of a 
sample. It will follow from the analysis below that the ba- 
sic properties of the CSS take their "canonical" form for 
the periodic boundary conditions (PBC). Other bound- 



ary conditions induce some characteristic features in the 
CSS which depend on the properties of the boundary. 

In order to study the qualitative role of the system ge- 
ometry for the CSS we consider the system which shows 
all the typical features of the critical system and allows 
for a rigorous analysis. This is a 2D weakly disordered 
electron system. The small parameter which allows for a 
rigorous treatment is the inverse dimensionless conduc- 
tance g^ 1 — A/E c <C 1, where h/E c — to = L 2 /D is 
the diffusion time and D is the diffusion coefficient. For 
not too large system size L <C £, where £ = ^exp[g] or 
£ = Zcxp[g 2 ] is the localization radius in the orthogonal 
or unitary ensemble, respectively, one can neglect the lo- 
calization effects and consider the dimensionless conduc- 
tance g as independent of the size L. At such conditions 
the behavior of the 2D system is similar to the one in the 
true critical point. 

The two-level correlation function R(s) in such 2D sys- 
tems has been considered in Ref. for the case of the pe- 
riodic boundary conditions only. It has been shown that 
the two-level correlation function R(s) is exponentially 
small at s > j if one neglects the weak-localization ef- 
fects. This means that the so called Altshuler-Shklovskii 
§| tail R(s) = C d s- 2+d / 2 is absent in 2D systems. This 
statement is true || for a generic critical system. 

In the present Letter we consider R(s) at s ^> g for 
an arbitrary 2D surface S of the area As- We will show 
that in the case where the surface has a non-transparent 
boundary, the two-level correlation function has a power- 
law tail R(s) ~ g _1 / 2 s -3 / 2 . Moreover, we will show that 
even in the absence of boundaries (e.g. for the sphere) the 
correlation function R(s) still has a power-law tail R(s) = 



C& 



This tail at s ^> g is of the same form as the 



universal WD tail for 1 -C s <§; g. However the coefficient 
C depends on topology of the surface and is equal to zero 
only for the torus topology which is equivalent to the 
PBC. 

The two-level correlation function is defined as: 



R(u) = \ 
P 
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where (...) denote averaging over all realizations, p = 
(p{E)) is the average one-electron density of states 
(DOS), and p{E) is the DOS for a particular realization 
of disorder: 



p{E) = S- 1 Y,S(E-E n ) 



(2) 
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where E n are exact one-electron energy levels. Since we 
are interested in the energy scale w = sA which vanishes 
in the TD limit, the average DOS is considered energy- 
independent. 

It is remarkable that in the limit g 3> 1 and s 3> 1 the 
spectral correlation function R(s) of the quantum prob- 
lem can be expressed [|| entirely in terms of the eigen- 
values e^A/h of the classical diffusion problem with the 
Neumann boundary conditions which correspond to the 
non-transparent boundary: 



R(s) 



1 
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where (3 = 1, 2, or 4 for the orthogonal, unitary and sym- 
plectic ensembles, respectively [||. Note that the dimen- 
sionless conductance g — e\ — eo is just the gap between 
the lowest (eo = 0) and the first non-zero eigenvalues. 

This means that independently of the details of the 
electron conduction band structure and the short-range 
correlated random potential, at g 3> 1 there exists a re- 
gion of s where the electron level correlations depend only 
on the spectrum of the Laplace-Beltrami operator Ac; 
on a curved surface: 



a c = — — (vgg ij — 

G ~ dx l dxi 



(4) 



where Q % i = [<j -1 ]ij, Q = detC?, and Q is a metric tensor 
on the surface. 

Eq. (||) can be rewritten in the following form ^| : 



K(t) 



' " ds _ Us „, , 1 
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e- Us R(s) = ^\t\p(\t\), (5) 



if one introduces the classical return probability P{t, r) 
to the point r: 



P(t,r) = 5^exp[-te M ] [$ M (r) 



(0) 



In the diffusion region the main contribution to the 
return probability is given by the continuous approxima- 
tion in Eq. (Q) . In this approximation p{t) is well known 
to be proportional to l/t in 2D systems thus giving rise 
to a constant K(t). It is the property of all the critical 
systems that the main contribution to Kit) is indepen- 
dent of t for sufficiently small times pXj ) pJ| . This means 
that the tail in R(s) for s 3> g is determined in critical 
systems entirely by the corrections to the return prob- 
ability. There are two types of such corrections. One 
source of corrections is the weak-localization effects con- 
sidered in Ref. or the finite-size corrections to scaling 
considered in Ref. p0|. In 2D systems they lead to: 



Rwl(s) 



T9 2 s \ 

-g- 3 s-\ 



0=1,4; 
[3 = 2 



(8) 



Another source of corrections which has not been con- 
sidered so far is the corrections to the continuous ap- 
proximation in Eq.(^) that is the difference between the 
sum and the corresponding integral. This correction is 
sensitive to the boundary conditions and topology of the 
surface. 

The sensitivity of pit) to the boundary conditions orig- 
inates from the r-dependence of the return probability 
P(t,r). If the point r is at a sufficiently small distance 
x <C L from a smooth boundary, the r-dependent cor- 
rection 8G[t, r, r) = SPb(t,r) to the Green's function 
of the diffusion equation with the Neumann boundary 
conditions, is given by the "image source": SP b (t,r) = 
(iirDT)- 1 exp[-a; 2 / 'DT}. Then assuming the smooth 
boundary of the length C one obtains ]12fl for the cor- 
responding correction to the averaged return probability 
p(t): 



Sp b (t)= / 5P b (t,r)Cdx = 



C 



(9) 



The boundary-induced tail in R(s) is found immediately 
from Eqs.(|),(|): 



and the averaged return probability: 



R b (s) = 



Pit) 



J s drP(t,r) =^>xp[- 



ts. 



(7) 



where ^ (r) is an eigenfunction of the diffusion (Laplace- 
Beltrami) operator and t = T/tjj = TA/h is time in 
units of the Heisenberg time tn ■ 

Eq.(||) sets the relationship between the tail in R(s) at 
s > 1 and the small-time behavior of the return prob- 
ability p(t) at t <C 1. There are two regions of t with 
different behavior of p(t) : the diffusion region t <C 1/g 
and the ergodic region l/j<f« 1. In the ergodic re- 
gion te M 3> 1 for all /i 7^ 0, so that only the zero mode 
£o = contributes to Eqs.(^),(^). In this region p(t) = 1 
is independent of t, and one obtains the universal WD 
result R(s) = —-p-p s~ 2 . 



8/3 (irgt) 1 / 2 {2irs) 3 / 2 



(10) 



where g b = D/C 2 A. Note that for £~Lwe have g b ~ g. 

Comparing Eqs.(p|),(|lC|) one can see that the 
boundary-induced power-law tail in R(s) is larger than 
the localization- induced tail if g <C s <C g 3 for /3 = 1,4 
and g <C s <C g 5 for (3 — 2. 

If the boundary is absent at all (R b (s) = 0), still 
there is a correction to the continuous approximation in 
Eq. (0) . This correction is universal and depends only on 
the topology of the surface. In order to understand how 
such a topological correction may arise at times T which 
are small compared with the diffusion time, let us con- 
sider the return probability P(t,r). For times T <C to 
(t <C l/<?) the diffusing electron probes only the vicinity 
of the starting point r. Therefore, the correction to the 
return probability P(t, r) may depend only on the local 
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curvature. We will show below (see Eq.©) that it is 

proportional to the Gauss curvature K r = k^ k^ in a 
point r, where k^ and k^ are the principal curvatures. 
An information about topology of a surface arises be- 
cause of the integration over all positions of the starting 
point r in Eq.(Q). Indeed, according to the Gauss-Bonnet 
theorem jl3| the integral of K r over a surface S is related 
with the Euler characteristic x(&) of the surface: 



/ K r dcr = 2nx(S). 
Js 



(11) 



The quantity x(S') = V + F — E is related with the num- 
ber of vertices V, edges E and faces F of the surface 
triangulation and depends on the connectivity of the sur- 
face, i.e. on the number p of "handles" (holes): 



X (S) = 2(l-p), 



(12) 



where p = for a sphere and p — 1 for a torus. 

In the simplest case of a sphere where each eigenvalue 
£fi = 1(1 + 1) of the Laplace-Beltrami operator Eq.(Q) 
is (21 + 1) times degenerate, the topological correction 
to p(t) at tg <C 1 can be obtained by a straightforward 
summation in Eq.(Q) using the well known |l4j| formula 

E~/(0*n /a /(0<fl + H/'(-V2): 



p(t) = J2(2l + 1) exp[-tgl(l 



1=0 



■K4 



(13) 



One can see that for a sphere (p — 0) the topological 
correction to p(t) at ij < 1 is a universal constant -|. 
Then using Eq.(|l2|) and <5p(i) oc x(S) we conclude that 
for a generic surface the topological correction reads: 



8ptop{t) 



1-p 



(14) 



This correction is time-independent and one immediately 
obtains from Eq.(|5|) a topological tail in R(s): 



Rtop(s) 



1-P 

3tt 2 /3s 2 



(s » g). 



(15) 



Comparing this correction with the weak-localization 
correction Eq.(||) we see that the latter is smaller at 
g < s < g 2 or g < s < g 3 for (3 = 1,4 or (3 = 2, 
respectively. Note that the topological tail has the same 
form as the universal WD tail. However, the former is 
valid in the diffusion region s 3> g while the latter is valid 
in the ergodic region s <C g. 

The result Eq.(|l4]) follows from the theory of the 
Laplacian on Riemannian manifolds (TLRM) |15| and 
has been used (l^] to find corrections to the semiclassi- 
cal density of states in quantum billiards and resonators 
p6[ . Here we present an elementary derivation of Eq. fll4| ) 
which does not require the full power of the TLRM. 

For an arbitrary non-singular point r on a surface, it is 
always possible to choose the local system of co-ordinates 



(Ci i €2) such that the inverse metric tensor G^i^i,^) m 
the vicinity of the origin (point r) takes the form: 



g^^g 22 ^i } g 12 = g 21 = -K r ^ 2l 



(16) 



where K r is the Gaussian curvature at a point r, and 
higher order terms in £ 1,2 are omitted. Then we obtain 
for the Laplace-Beltrami operator Eq.fl): 



d 2 



d 2 



di 2 d% 

= A + V. 



d d d d 

Wi W2 + W2 Wi 



(17) 



As has been mentioned above, for small times T <C to 
the term proportional to the Gaussian curvature K r 
can be treated as perturbation. The return probabil- 
ity P(t, r) is given by the Green's function G(t,T, r) 
of the diffusion operator d/dt — Ag, where t = DT. 
Then all we should do in order to find the correction 
to the local return probability is to compute the correc- 
tion to the Green's function G (1) = G^VG {t) \ where 
G( 0) = (471-t)- 1 exp[-(£ - S,') 2 /At]. The result is: 



P(t,r) 



As_ K*_ 

AttI 12tt' 



(18) 



where ^5 is the surface area. The first term in this ex- 
pression corresponds to the continuous approximation. 
The second one is the correction due to the local cur- 
vature K r . This very term has been used in deriving 
Eq.©. 

In conclusion, we have considered the power-law tail in 
the two-level correlation function R(s) for 2D disordered 
systems with the diffusion motion of electrons. In the 
limit of large dimensionless conductance g 1, there is 
an interval of s 3> g where the tail is entirely determined 
by the geometry of the sample. It consists of the two 
contributions: a boundary contribution given by Eq.([lO|), 
and a topological contribution given by Eq. (|l5| ) which are 
the main results of the Letter. For closed surfaces with- 
out a boundary, the tail in R(s) for s 3> g is determined 
only by topology and is of the same form as the uni- 
versal Wigner-Dyson tail for 1 « s « j. However, the 
numerical topological pre-factor depends on the connec- 
tivity of the surface and is negative for single-connected 
surfaces (sphere), zero for surfaces of the torus topology, 
and positive for surfaces of higher genus. 
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